Protein crystals and charged surfaces: interactions and heterogeneous nucleation 
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As proteins typically have charges of around 10, they will interact strongly with charged surfaces. 
We calculate the electrostatic contribution to the interaction of crystals of protein with charged 
surfaces. The surfaces repel like-charged crystals and attract oppositely-charged crystals, with free 
energies which can be easily several kT per protein molecule brought into contact with the surface. 
This means that oppositely charged surfaces can act as a nucleant, they can induce nucleation of a 
protein crystal by lowering the free energy barrier to heterogeneous nucleation of the crystal from a 
dilute solution. 



I. INTRODUCTION 

Here we will consider the interaction between charged 
surfaces and protein crystals. Proteins are themselves 
charged and so we would expect them to interact strongly 
with a surface that is charged. Using the Poisson- 
Boltzmann equation we find that indeed protein crystals 
are attracted by oppositely charged surfaces and repelled 
by surfaces with charges of the same sign. In itself this is 
hardly surprising but with quite simple calculations we 
quantify this attraction and repulsion. We show that eas- 
ily accessible surface charge densities of the opposite sign 
to that of the protein molecules, are able to greatly favour 
heterogeneous nucleation at the surface by reducing the 
free energy of a nucleus at the surface by several faT 
per protein molecule at the surface. Controlled hetero- 
geneous nucleation is vital to the production of protein 
crystals, which are required for structure determination 
via X-ray crystallography. 

At a first-order phase transition, such as the crystalli- 
sation of a protein from solution, the transition starts 
with nucleation [U0. Protein crystallisation is the 
"main bottleneck" Q in the determination of the three- 
dimensional structure of proteins. Determining this 
structure is crucial for understanding what a protein docs 
and how it does it. Protein crystallographers wish to: i) 
induce nucleation in the relatively weakly supersaturated 
solutions, within which protein crystals grow slowly and 
so incorporate few defects, and ii) control heterogeneous 
nucleation, in particular have nucleation occurring only 
at specific locations on a surface. If crystals only nucleate 
on widely separated patches on a surface then the growth 
of one crystal will not interfere with and limit the growth 
of another crystal. 

See for example Refs. 0, Q for experimental work on 
adding solid surfaces to protein solutions to induce nu- 
cleation, and see Refs. ja, Q for experimental work on 
patterned surfaces, which shows that protein crystals can 
be made to grow preferentially on specified patches of a 
surface. In Refs. H a patterned surface is used in 
which patches of the surface are one form of doped Si 
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and the rest is a different form. The protein nucleates 
preferentially on one of these forms and so by controlling 
where this form is found in the pattern, the crystals can 
be made to nucleate and grow far from each other, fa- 
cilitating the formation of large crystals. Here we show 
that a pattern of two different surface types of opposite 
charges can produce sufficiently large differences in sur- 
face free energy to make heterogeneous nucleation much 
easier on one surface than on the other. 

We note that here we only consider the effects of charge 
interactions, screened by a 1:1 electrolyte. Other inter- 
actions, such as short-range specific interactions between 
the protein molecules and surface are known to be impor- 
tant, see the review of Ostuni et al. 8]. When these in- 
teractions are important, our calculations, which neglect 
them, will only yield estimates of the variation with 1:1 
salt concentration, of the interaction free energy of a pro- 
tein crystal or crystalline nucleus with a surface. How- 
ever, experimental systems do exist with high, variable 
charges and highly hydrophilic surfaces. These highly 
hydrophilic surfaces will have only weak non-electrostatic 
interactions with the proteins, they also has the advan- 
tage of minimising the problem of protein denaturation. 

One example is a highly charged self-assembled mono- 
layer, such as has been shown to adsorb polyelectrolytes 
p|, another is a membrane containing charged lipids 
[ToL ITU . Both of these systems should be well described 
by the theory presented here. There is a great deal of 
work on the association of proteins and DNA with op- 
positely charged lipid membranes, Refs. 0, |n| are just 
two examples of a large literature. As far as the author is 
aware, there has been no experimental study of adsorp- 
tion of protein on a highly charged, highly hydrophilic, 
self- assembled monolayer, but such a system would offer 
the opportunity to study the effect of electrostatics with 
minimal interference from other interactions. This is true 
both for the adsorption of single proteins and of protein 
crystals. 

Our calculations will all be for infinite planar inter- 
faces, despite the fact that our primary interest is in 
the nuclei of protein crystals against charged surfaces, 
not bulk crystals against surfaces. Although, nuclei are 
small, 10s of protein molecules, the part of the free energy 
which scales with the contact area between the nucleus 
and the surface will still dominate the edge contribution, 
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FIG. 1: A schematic of a nucleus of a protein crystal against 
a charged surface. The hatched circles are the proteins in the 
crystal lattice, each with its charge Q marked. The surface is 
shaded grey and the surface charges, positive in this case, are 
shown as +'s. 
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as the Debye length is typically a few nms at most. Thus, 
we consider only the dominant surface term, consider- 
ing sub-dominant terms within simple theories of simple 
generic models is not useful. We will also leave consider- 
ation of defects in the surface to later work. 

In the next section we calculate the electrostatic con- 
tributions to the surface free energy. The third section 
includes example results and a discussion and the last 
section is a conclusion. 



II. CALCULATION 

Here we calculate the electrical double layer contribu- 
tions to the interfacial free energies of solid-salt-solution, 
solid-protein-crystal and solution-protein-crystal inter- 
faces. The calculation for the solution-protein-crystal 
interface is just a repeat of that done in Ref. ^3] and so 
will not be described in detail. We will use the model of 
Ref. 01 f° r the protein crystal and solution, along with 
a simple model of a charged solid surface. As we are in- 
terested in interfacial free energies it is most convenient 
to work at constant chemical potential of the salt, and 
so the appropriate free energy is the grand potential f2. 
We will split this grand potential into the part due to 
the charges on the surface and on the protein molecules, 
together with the salt and counter ions, Q e i, and another 
part which includes the rest of the interactions, O ne z, 

il = Quel + &el- (1) 

fl ne i includes effects such as the dispersion interactions 
between the protein molecules and the surface, and any 
short-range interactions between the surfaces of the pro- 
teins and the solid surface. We will not calculate fi ne ;, 
and so will only be able to estimate the absolute values of 
the interfacial free energies, when n ne i is much smaller 
than fl e i, for example with highly charged, highly hy- 
drophilic, surfaces, and highly charged proteins. How- 
ever, if we assume that fl ne i is weakly dependent on the 



salt concentration, and that the charges on the surface 
can be varied without varying Q ne i significantly, our re- 
sults will describe the variation in the interfacial free en- 
ergies with surface charge density and the salt concen- 
tration. 

Rather than considering a specific protein crystal, we 
take over the jellium concept from the theory of metals 
to make a general estimate of the effect of salt on the in- 
terfacial tension. We replace the detailed charge density 
due to the protein molecules by a uniform background 
char ge d ensity p p , cut off abruptly at the interface. In 
Ref. j Sear and Warren did so for the protein-crystal- 
salt-solution interface, and this followed on from the work 
of Warren ^| on the bulk phase behaviour. Here when 
we calculate the interfacial free energy for a protein crys- 
tal at the surface, we will assume the protein crystal is 
flush against the surface and so the charge density is uni- 
form right up to the surface. Our model of the charged 
surface of the substrate is that of an infinite charged plane 
at a fixed charge per unit area a. The solution outside 
the crystal is taken to be sufficiently dilute that we can 
neglect any protein molecules in it and treat it as a salt 
solution at a concentration p s . For defmiteness we take 
the protein to be positively charged so the counterions 
are negative and the coions are positive; both are mono- 
valent. The surface is in the xy plane at z = 0, with the 
solid in the z < half-space and either a salt solution or 
a protein crystal in the z > half-space. 

We will use the same notation as in Ref. [l^, so we 
work in units where e = k^T = 1. In these units, the 
Coulomb potential energy U between a pair of elemen- 
tary charges separated by r is U = ls/r , where Ib is the 
Bjerrum length, equal to 0.72 nm in water at room tem- 
perature (7b = e 2 /47refceT). We assume a constant value 
of Ib, and ignore dielectric effects. The counterion and 
coion densities are p-[z) and p+(z), respectively. In the 
bulk salt solution, both tend to the salt concentration p s , 
and inside the protein crystal (p_ — p + ) — » p p . 

We will use a grand potential, f2 e / , which contains only 
ideal solution terms for the ions and the associated elec- 
tric field at the interface. The grand potential f2 e /, is 

m, 

roc p2 

Q el = dzcj(z), 0; = V] ft(ln — - 1) + . ( 2 ) 
Jo Ps ottIb 

The first terms in u> are the ideal solution terms (the ions 
share a common chemical potential p± = lnp s ). The last 
term is the electrostatic energy, wherein E = —d<p/dz is 
the electric field strength corresponding to an electro- 
static potential cf> which satisfies the Poisson equation, 

d 2 (j> _ ( -47rZ B (p+ — p~) salt solution 

dz 2 ~ \ -47T? B (ft+ - P- + P P ) protein crystal ' 

The variational principle SQ e i/8[p±(z)] = applied in 
this problem yields p±(z) — p s exp[=p</>]. The elec- 
trostatic potential then satisfies the Poisson-Boltzmann 
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equation, 



dz 2 
where 



k c sinh 6 = 



salt solution 

—AttIbPp protein crystal 



«2 



8ttI b Ps 



(4) 



(5) 



For our assumption of a planar surface in the z = xy 
plane, with a fixed charge per unit area <r, the boundary 
condition on at z = 0, is that the electric field at the 
surface 



4ttIb<7. 



(G) 



The other boundary condition is that d<fi/dz — 0, z — > oo 
A. Surface— salt-solution interface 



If Eq. is not linearised the potential is given by 
El = 21n[(l + C*e" K ^)/(l - Ce- K ° z )} where C 
is determined by the boundary condition Eq. JSJ and is 
given by 



C 



2ttIb<J 



(10) 



Thus, we can evaluate the surface free energy once we 
have converted Eq. JSJ into an expression in terms of 
only. Equation (0J in the salt solution can be integrated 
once with respect to 0, 



\ dz 



I — I = 2/^(cosh0-l) 



(11) 



which gives us the electric field term in to. The other, 
ideal gas, term can be written in terms of p s and using 
P±( z ) = Ps exp[=F0]. Then A7 S is given by the integral 



This is just the problem addressed by Gouy and Chap- 
man at the beginning of the last century 14]. For sim- 
plicity we will only derive the expressions in the regime 
where the equations may be linearised. Then we will sim- 
ply state the general result obtained when the equations 
are not linearised. Linearising the Poisson-Boltzmann 
equation, Eq. (@J, inside the salt solution, and using the 
boundary conditions, we have 



(z) = 4nlBcrK s 1 exp(— n s z) z > 0. 



(7) 



A simple exponential decay from a value at the surface, 
proportional to the surface charge. 

The surface free energy is, by definition, the difference 
between the actual grand potential per unit area of the 
surface-salt-solution interface and that it would have if 
the salt solution continued unperturbed right up to solid 
surface. We therefore have to calculate the grand poten- 
tial, Eq. (J2J, then subtract the grand potential for the 
bulk solution: Eq. with the actual ui(z) replaced by 
its value in the bulk salt solution, oj(+oo) |15|| . 



A7s 



dz (to(z) — aj(+oo)) . 



(8) 



This is a general expression, in the linear regime we can 
expand out u> of Eq. @ and keeping only the terms up to 
quadratic order in 0, we have that lu = — 2p s +2p s <f) 2 . The 
bulk value cj(+oo) = —2p s , as the potential in the bulk 
of the salt solution is taken to be zero. The surface free 
energy is then obtained by substituting these expressions 
into Eq. (jHJ. Doing this and using Eq. Q), we have 



A7 S = 2ttIbc k 



2-1 



linear, 



(9) 



which is positive and proportional to the square of the 
surface charge and to the screening length, as we might 
have expected. 



A7 S = 2p s / G?z0sinh0, (12) 



which is straightforward to evaluate numerically. 



B. Surface— protein-crystal interface 

As with the surface-salt-solution interface we start 
with the linear Poisson-Boltzmann equation. Here the 
potential is that inside the protein crystal so we linearise 
Eq. (0J and obtain 



f0 
dz^ 



Pp 
2 Ps 



(13) 



where we used Eq. © for n s . The second term inside the 
parentheses is minus the Donnan potential 0d inside the 
protein crystal. The general expression for the Donnan 
potential is 



sinh < 



P P 
2 Ps ' 



(14) 



which when linearised gives 0d = Pp/{2p s ). Thus, the 
solution to Eq. Q13f) . with the boundary condition of fixed 
charge density, Eq. ©, is 



(z) = 0d + 47rZBcr/t s 1 exp( 



z > 0. 



(15) 



In the bulk protein crystal, = 0d ■ We expand the u) 
of Eq. in powers of 0d and keep only the terms up to 
quadratic order, and find that in the crystal u>(— oo) = 
— 2p s + 2p s 0|, 12]. Performing the same expansion at 
the interface, we obtain for the excess grand potential at 
the interface 



w(-oo) = p s 



(16) 
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where we used Eq. (|15|l to obtain the derivative d(j>/dz. 
The interfacial free energy per unit area A7 X is the inte- 
gral over all z of the excess grand potential of Eq. iJTfljl . 
It is 



A 7x = 2ttIbc 2 k s 1 + (7 0rj linear, 



(17) 



If we do not linearise the Poisson-Boltzmann equation, 
Eq. we can still integrate it once to obtain 



FIG. 2: The electrical double layer contribution to the change 
in free energy per unit area, times the surface area of the 
crystal per protein molecule, when a protein crystal is brought 
from the bulk of the salt solution to the charged surface of a 
solid, A/ 2 , Eq. 1231 . The solid, dotted, dot-dashed and dashed 
curves are for charge densities ae = 1, 0.1, —0.1 and — lnm -2 , 
respectively. Results are shown for salt concentrations down 
to 0.025M. In the limit of the salt concentration tending to 
zero the surface free energies diverge. 



^ I = 2n 2 (cosh(/> 
dz ' 



cosIi^d) — 87tZbPp i 



<h)> 

(18) 

Combining this with the boundary condition, Eq. © we 
obtain an equation for the potential at the surface <p(z = 
0), 

[2k 2 (cosh. 4>(z = 0) — cosh^n) — 8ttIbP p (<p(z = 0) — <frn)] 



1/2 




= ±47tZbCT, 

which can be solved for <j)(z — 0), and then once this is 
known the profile 4>(z) is readily obtained by numerically 
integrating Eq. (|TH|l . On the right hand side the + (— ) 
sign is taken when a > (a < 0). The interfacial free 
energy is then obtained from 



A 7x 



dz (lo(z) — u>(— oo)) , 



(20) 



where the grand potential per unit volume in the protein 
crystal is u>(— oo) = — 2 j o s (cosh </>d — (f>u sinh (f)B>) ■ Using 
p±(z) = p s exp[=F</>], and Eq. ljl"B|) for the electric field 
term in Eq. @, we have that the electrical double layer 
contribution to the free energy of the solid-protein inter- 
face is given by 



A 7x = 2p s I dz 0[sinh0 — sinhi 
Jo 



(21) 



C. Salt-solution— protein-crystal interface 

This interface was the subject of Ref. I* 1 the high 
salt regime where the Poisson-Boltzmann equation can 
be linearised the electrical double layer contribution to 
the free energy of the interface between a protein crystal 
and a salt solution is 



A 7e l 



Ps4>T> 



linear. 



(22) 



See Ref. 121 for the full nonlinear calculation. 



III. RESULTS 

We can calculate the free energies of the three inter- 
faces: solid-salt-solution, solid-protein-crystal and salt- 
solution-protein-crystal. This allows us to see whether 



(19) 
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a protein crystallite, such as a nucleus or growing crys- 
tallite in solution, is attracted or repelled by the surface. 
It is attracted if the free energy change on bringing a 
protein crystal from the bulk solution to the surface is 
negative and repulsive if the change is positive. When 
a protein crystallite is brought from the bulk of the so- 
lution into contact with the surface two interfaces are 
destroyed, that between the salt solution and the protein 
crystal and that between the salt solution and the solid, 
and one is created, that between the solid and the pro- 
tein crystal. Thus, the free energy change per unit area 
of protein crystal brought into contact with the solid is 



A = -A 7s - A 7cl + A 7x , 



(23) 



which in the high salt limit is, using Eqs. ©, l|17|) and 

(23 > 



A = 



Ps<t>B 

2k« 



linear. 



(24) 



The first term comes from the disappearance of the salt- 
solution-protein-crystal interface and is positive as the 
electrical double layer contribution to the interfacial ten- 
sion of this interface is always negative; see Ref. [l2j for 
a discussion. As <^>d has the same sign as that of the 
charges on the protein molecules (here taken to be posi- 
tive) , the second term is positive if the surface and pro- 
tein molecules have charges of the same sign and negative 
if the have opposite signs. As we should expect, surfaces 
repel crystals of proteins with the same charge. They 
attract crystals of oppositely charged protein molecules, 
providing that the surface charge a is sufficiently high. 
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If the charge density is small or zero then they repel pro- 
tein crystals due to the free energy cost of destroying the 
salt-solution-protein-crystal interface. 

Equation (|24[1 is obtained by linearising the Poisson- 
Boltzmann equation. We have also solved the full equa- 
tion and present example results in Figs. [21 and [21 The 
parameters for the protein are those used previously |l2j | , 
and were chosen to model lysozyme. This is a small pro- 
tein with quite a large net charge for its size, whose be- 
haviour in NaCl solutions has been extensively studied 
[lit ITfil IT7L ITsl IT^ | . The model used here has been to 
shown to give the variation with NaCl concentration of 
the solubility of lysozyme crystals correctly to within a 
factor of about 2, it overpredicts the variation |l3l 
The agreement between its predictions and the surface 
force apparatus measurements of Sivasankar et al. [2(j 
of the repulsion between monolayers of the protein strep- 
tavidin, is comparable Jl 2 | . To model lysozyme we took 
a charge Q = 10 0, Il7ll2l| , and a charge density due to 
the protein p p = 0.25nm -3 , obtained by taking the vol- 
ume per lysozyme molecule in the crystal to be 40nm 3 . 
The surface area of the protein crystal per lysozyme 
molecule is taken to be I 2 = 12nm 2 [T^j. For proteins 
such as lysozyme which have a significant net charge we 
would expect our predictions to be approximately as ac- 
curate as those of earlier predictions of this theory, i.e., 
predicting the trends correctly but giving numbers which 
are a factor of 2 or more out. We expect our predictions 
to be unreliable if the net charge is small, then the inter- 
action between the dipole moment of the protein and the 
surface may dominate, or if the charge is distributed on 
the surface in a way which is highly inhomogeneous. If 
the charge density within the volume within a few Debye 
lengths of the surface differs greatly from p p , then our 
assumption of uniform charge density due to the protein 
will be poor. See Ref. [12| for further discussion of the 
assumption that the protein crystal can be modeled by a 
step-function charge density. 

In Fig. [21 we see that for lysozyme and solid surfaces 
with a = ±1 nm -2 , the free energy change of bringing 
the protein crystal into contact with the solid is substan- 
tial unless the salt concentration is around 1M. Even 
at a salt concentration of 1M, the free energy change is 
close to 2ksT per lysozyme molecule at the surface of 
the protein in contact with the solid surface. Whereas 
at salt concentrations of about 0.1M and lower the free 
energy change is of order 10/cbT per lysozyme molecule. 
Note however, that we have assumed that the protein 
crystal is in a salt solution, i.e., a solution with negligible 
amounts of protein in it. Decreasing the salt concentra- 
tion increases the solubility of lysozyme , and so this 
assumption will tend to worsen as the salt concentration 
drops [12L [l3| . The salt is not only changing the inter- 
action between the protein and the surface but also the 
interaction in between protein molecules and this needs 
to be borne in mind. 

In Fig. [21 we see that at high salt concentrations, 1M, 
the free energy change on bringing a crystal into con- 



FIG. 3: The electrical double layer contribution to the change 
in free energy per unit area, times the surface area of the 
crystal per protein molecule, when a protein crystal is brought 
from the bulk of the salt solution to the charged surface of a 
solid, Ai 2 , Eq. 1231 . The solid, dashed and dotted curves are 
for salt concentrations of p s = 0.1, 0.5 and 1M, respectively. 
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tact with the surface is never more than a couple of k^T, 
even for quite highly charged surfaces. However, for more 
modest salt concentrations of 0.1M, rather large free en- 
ergy decreases per protein molecule are easily obtained at 
modest charges per unit area on the surface. Also, note 
that A is slightly positive for an uncharged surface. This 
is due to the fact that the electrical double layer con- 
tribution to the protein-crystal-salt-solution interface is 
negative, Eq. I)22[l. and this contribution is lost when the 
protein crystal is brought into contact with the surface. 

The potential near the surface, both when the salt so- 
lution is in contact with the surface and when the protein 
crystal is, is shown in Fig.^ When the protein crystal is 
in contact with the surface the potential tends to its value 
in the bulk of the crystal: the Donnan potential, which 
is 1.48 here. The potential at the surface is larger than 
one, and so the linear approximation is quite poor for the 
charge a = ±1 nm~ 2 . In particular, Eq. (|24|l significantly 
overestimates the magnitude of A at high surface charges. 
However, Eq. I|24(l is reasonably accurate for the smaller 
charge per unit area, a = ±0.1 nm -2 . As an example, 
for a = -lnm~ 2 and p s = 0.1M, Al 2 = -11.1 without 
linearisation, whereas Eq. 1|24|) predicts AZ 2 = —17.0. 



IV. CONCLUSION 

We have used an existing [T2I generic model of 
a protein crystal to estimate the electrostatic contribu- 
tion A to the free energy change when a protein crys- 
tal is brought into contact with a surface with a fixed 
charge density a. Unsurprisingly, the free energy change 
is negative for a surface with a charge opposite in sign 
to that on the protein molecule, and positive if the sur- 
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FIG. 4: The potential, in units of k^T/e, as a function of 
distance z from the surface. The solid curves are for a salt 
solution in contact with the surface and the dashed curves 
are for a protein crystal against the surface. The upper solid 
and dashed curves are for a positively charged surface, a = 
lnm -2 , and the lower curves are for a negatively charged 
surface, a = — lnm -2 . The salt concentration is 0.1M, and 
the protein is positively charged, p v — 0.25 nm -3 . 
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face has a charge with the same sign. We quantified 
this and found that very reasonable charge densities, of 
1 elementary charge per nm 2 or less, were sufficient to 
achieve free energy changes per protein molecule at the 
surface of several k^T or more, provided the salt con- 
centration was less than around 1M. The electrostatic 
contribution to the free energy is the energy associated 
with the electrostatic interactions in between the sur- 
face, protein molecules, and salt and counterions, and the 
translational entropy of the salt and counterions. Typi- 
cally, there will also be non-electrostatic interactions of 
the proteins with the surface, particularly if the surface is 
at least moderately hydrophobic. Then the proteins may 
even unfold resulting in irreversible adsorption. Thus, 
our results will be relevant only to those surfaces which 
are sufficiently hydrophilic that the proteins interact only 
weakly with the surface and do not unfold. However, ob- 
taining such surfaces is relatively straightforward via the 
use of self-assembled monolayers 0, 0. Our model of 
the protein is a very simple one, appropriate for proteins 
with reasonably large net charges, so that the terms we 
calculate should be dominant, and where this charge is 
distributed so that the charge density facing the surface 
is neither much higher nor much lower than the average 
charge density on the surface of the protein. 

Perhaps the most important reason for considering 
protein crystals at surfaces is interest in their hetero- 
geneous nucleation. Protein crystals nucleate at surfaces 
and so without an understanding of crystallites at sur- 
faces we cannot hope to understand how the crystals ac- 
tually form. Crystallising proteins is essential for the 
determination of their all-important three-dimensional 



structure via X-ray crystallography Heterogeneous 
nucleation is an activated process |j| and so proceeds at 
a rate which scales as exp(— AF* /k^T), where AF* is 
the height of the free energy barrier which must be over- 
come. Typical sizes of the critical nuclei in nucleation 
are a few 10s of protein molecules, and so one face of a 
crystalline nucleus has perhaps 10, or a few less, protein 
molecules. The critical nucleus is that at the top of the 
barrier, the nucleus which requires a free energy A_F* to 
create it pj. Thus, the variation in AF* with A, will be 
approximately 10 Al 2 . From Figs. [21 and |3 we see that 
this may easily be 10s of k^T, causing increases or de- 
creases in the rate of heterogeneous nucleation of many 
orders of magnitude. The barrier to nucleation may be 
lowered by 10s of ksT, inducing nucleation in solutions 
which would otherwise be metastable, or it may be raised 
by 10s of k^T preventing nucleation from occurring on 
one part of the surface. The former is required if a sur- 
face is to be used as a nucleant: a surface which triggers 
nucleation 0,0. The latter is required for spatial control 
of heterogeneous nucleation 

Finally, we consider wetting by the crystalline phase 
of the surface-salt-solution interface. At coexistence, for 
a bulk solution phase in contact with the surface, a di- 
rect surface-salt-solution interface is not the only pos- 
sibility. It is possible for a slab of the protein crystal 
to be interposed between the surface and the solution 
phase. The direct surface-salt-solution interface is then 
replaced by by a surface-protein-crystal interface plus 
solution-protein-crystal interface. The protein crystal is 
said to have wet the surface-salt-solution interface. Wet- 
ting occurs at equilibrium whenever it lowers the free 
energy 0, l^]. Denoting the full interfacial tensions 
of the surface-salt-solution, surface-protein-crystal and 
solution-protein-crystal interfaces by j s , "f x and 7, re- 
spectively, we have that the free energy change on inter- 
posing a slab of protein crystal between the surface and 
the solution is — "f s +Jx+J- Note that here the interfacial 
tension of the solution-protein-crystal interface appears 
with a positive sign whereas it appears with a negative 
sign in the definition of A, Eq. I|23[) . When a layer of the 
protein crystal is created so is a solution-protein-crystal 
interface, whereas when an existing protein crystal sur- 
rounded by solution is brought into contact with the sur- 
face a solution-protein-crystal interface is destroyed. 

Determining the sign of — 7 S +7^+7 is not possi- 
ble without determining both the electrostatic and non- 
electrostatic contributions to all three interfacial ten- 
sions. However, in Ref. |l2j| . the interfacial tension 7 
was estimated to be a few k^T per area of a protein 
molecule, i.e., jl 2 = 0(1). If we further assume that 
— 7s + lx is dominated by electrostatics, then in the lin- 
ear regime — 7 S + 7^ ~ c0d, from Eqs. @ and lfT7Jl . 
Thus for the electrostatic attraction between an oppo- 
sitely charged surface and a protein crystal to drive wet- 
ting of the surface-salt-solution interface, we require that 
|ct0d|^ 2 > 0(1). For highly charged surfaces, a of order 
lnm -2 , and low salt concentrations, this inequality is 
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easily satisfied. 

If the protein crystal wets the surface-salt-solution in- 
terface, then assuming that the layer of protein crystal 
can itself nucleate on the surface it will do so and the 
bulk protein crystal can grow from this layer at the sur- 
face. Then the barrier to crystallisation will be close to 
zero and the protein will readily crystallise from solution 
[23| . This assumes that the layer itself can nucleate, see 
Ref. |22| for an introduction to wetting layers, both at 



equilibrium and their nucleation. 
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